Quantum Trajectory Approach to the Stochastic Thermodynamics of a Forced 
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I formulate a quantum stochastic thermodynamics for the quantum trajectories of a continuously- 
monitored forced harmonic oscillator coupled to a thermal reservoir. Consistent trajectory- 
dependent definitions are introduced for work, heat, and entropy, through engineering the thermal 
reservoir from a sequence of two-level systems. Within this formalism the connection between ir- 
reversibility and entropy production is analyzed and confirmed by proving a detailed fluctuation 
theorem for quantum trajectories. Finally, possible experimental verifications are discussed. 



I. INTRODUCTION 

Thermal fluctuations cause the motion of a small clas- 
sical system, like a colloidal particle immersed in a vis- 
cous fluid, to be random and erratic. Such systems are 
clearly not within the scope of macroscopic thermody- 
namics Nevertheless, thermodynamic quantities - 
such as work, heat, and entropy - can be deflned consis- 
tently along individual stochastic trajectories within the 
theoretical framework of stochastic thermodynamics ^2- 
0] ■ Its two primary components are a flrst-law-like energy 
balance equation, introduced by Sekimoto in the context 
of stochastic energetics , and a definition of entropy 

along single fluctuating trajectories. The predictions of 
stochastic thermodynamics, which have been verified ex- 
perimentally |lll4l3| . have been important in categoriz- 
ing the fiuctuation theorems (reviewed in Refs. (6l.[l#IT7|) 
and in sharpening our understanding of thermodynamics 
at the nanoscale, especially with respect to the second 
law of thermodynamics [1, • 

Despite its significance, stochastic thermodynamics 
has yet to be extended to quantum mechanical systems. 
Like its classical counterpart, a quantum stochastic ther- 
modynamics would be beneficial for analyzing the fluc- 
tuations of thermodynamic quantities. In particular, a 
quantum stochastic thermodynamics could provide in- 
sight into quantum extensions of the work fluctuation re- 
lations jl^, [3 ■ ^ collection of predictions regarding the 
fluctuations in the work performed on a system driven 
away from equilibrium [T^, [13, [20l - [27j . Derivations of 
quantum work fluctuation relations begin with a defini- 
tion of work and a method to measure its fluctuations. 
For closed quantum systems - systems that do not ex- 
change energy with their surroundings [28| - the work 
performed by a quantum system during a thermody- 
namic process is determined by measuring the system's 
initial energy U and final energy U'; the difference is iden- 
tified as the work W ^ AU ^ U' - U \18., M'M- Stud- 
ies of various models have confirmed this quantum work 
fluctuation relation for closed quantum systems (36l - l40j 
and have pointed to its limitations 41j. However, these 
predictions have never been verifled experimentally, and 
only two experiments have been proposed [l^ . When 



the quantum system is open - exchanges energy with a 
thermal reservoir - there are a number of proposals for 
the definition of work. A common procedure is to mea- 
sure the energy of the thermal reservoir at the begin- 
ning e and at the end e' of a thermodynamic process in 
addition to measuring the system's initial and final ener- 
gies [l8l.[l9l. l43l - l45j . The change in the energy of the ther- 
mal reservoir is identified as minus the heat absorbed by 
the system, Q = — Ae = — (e'— e), and the work is inferred 
from the relation W = AU — Q. However, the possibil- 
ity of experimentally realizing a measurement of the en- 
ergy of an infinite thermal reservoir is still an open ques- 
tion [lH, . Esposito and Mukamel derived a quantum 
work fluctuation relation for an open quantum system 
defining the work along trajectories constructed from for- 
mal solutions of a quantum master equation [47] , though 
the connection between the work defined in Ref. [47[ 
and a measurable physical quantity remains unclear. In 
the high-friction limit, Deffner and Lutz investigated the 
work performed by an open quantum system [48-5C|. 
However, in this limit the equations of motion are effec- 
tively classical, as there are no quantum coherences (yet 
quantum fluctuations remain in addition to the thermal 
fluctuations). A quantum work fluctuation theorem has 
also been predicted by De Roeck and Maes for a path- 
dependent work defined along unitary quantum evolu- 
tions interrupted by projective measurements [51]. Yet, 
another approach was put forward recently by Hu and 
Subi§i using a decoherent histories aiialysis of a non- 
markovian quantum brownian motion [52']. Despite the 
many studies of work relations in open quantum systems, 
there are still unresolved questions regarding the possi- 
bility of experimentally measuring the work in a coherent 
quantum system, such questions that could be resolved 
by constructing a quantum stochastic thermodynamics. 

In this article, I take a first step in formulating a 
quantum stochastic thermodynamics using the quan- 
tum trajectory formalism, originally developed in the 
field of quantum optics in response to experiments on 
continuously- monitored individual quantum systems [s^- 
[5^ . Within this formalism quantum trajectories are the 
stochastic evolution of a quantum system conditioned on 
the outcomes of a sequence of weak measurements [s^ ]. 
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For concreteness, I analyze a specific model inspired by 
current experiments in cavity quantum electrodynam- 
ics [60l - [6^ : a forced quantum harmonic oscillator weakly 
coupled to a thermal reservoir. However, instead of con- 
sidering an infinite thermal reservoir, the thermal reser- 
voir is engineered by coupling the harmonic oscillator 
with a rapid succession of two-level systems, one at a 
time 63, 64]. By measuring the state of each two- level 
system after its interaction with the harmonic oscilla- 
tor, we are able to continuously monitor the evolution 
of the harmonic oscillator and measure the amount of 
energy transferred to the thermal reservoir, a strategy 
originally proposed by Crooks [i^ Derezihski, De 
Roeck, and Maes have also proposed the use of quantum 
trajectories in the study of quantum fiuctuation theorems 
for currents, but did not comment on their experimen- 
tal realizability [5l|, H^. By studying a concrete model, 
I am able to devise experimentally verifiable definitions 
of thermodynamic quantities along individual quantum 
trajectories. 

We will see that because the forced quantum harmonic 
oscillator's evolution is perturbed by its interaction with 
the sequence of two-level systems, its effective equation 
of motion is a stochastic Schrodinger equation. The 
derivation presented in this article will follow the ap- 
proach developed in Refs. [H, 113, |63, HI] , though adapted 
for a time-dependent Hamiltonian. This approach to 
open quantum systems, utilizing repeated interactions, 
has been studied rigorously by Attal and Pautrat [6^ 
and Attal and Joye |70| employing the quantum stochas- 
tic calculus devised by Parsatharathy [7l|. While the 
present analysis does not utilize Parsatharathy's quan- 
tum stochastic calculus, Pellegrini ^7^] and later At- 
tal and Pellegrini [t^ demonstrated that the stochastic 
Schrodinger equation considered here can be naturally in- 
corporated into this framework. Stochastic Schrodinger 
equations have also appeared in dynamic reduction mod- 
els, where they are postulated as a means to dynami- 
cally localize a state vector in Hilbert space [TJ]. The 
novelty of the present endeavor is to consider explicitly 
time-dependent Hamiltonians. By applying the secular 
approximation to the full time-dependent Hamiltonian 
(without utilizing Floquet theory as in Ref. ^595), I de- 
rive a stochastic Schrodinger equation and a quantum 
master equation valid even for nonadiabatic driving; a 
result that is necessary as we are interested far-from- 
equilibrium quantum thermodynamics. 

The construction of a quantum stochastic thermody- 
namics begins in Sec. HI] by specifying the model and in- 
troducing the stochastic Schrodinger equation. Then in 
Sec. lHIl I formulate a quantum stochastic thermodynam- 
ics by defining heat, work, and entropy along individual 
quantum trajectories. Work fiuctuation relations are an- 
alyzed in Sec. IIVI where a detailed fluctuation relation 
for quantum trajectories is derived. Finally, I conclude 
in Sec.|V]with a discussion of possible experimental verifi- 
cations of quantum stochastic thermodynamics in cavity 
quantum electrodynamics experiments [60l - l62| . 



II. FORCED QUANTUM HARMONIC 
OSCILLATOR COUPLED TO A THERMAL 
RESERVOIR 

In this section, I specify the details of the model and 
fix notation. The model is a one-dimensional forced har- 
monic oscillator weakly coupled to a thermal reservoir. 
Since a common experimental realization of the quantum 
harmonic oscillator is a single mode of an electroma gen - 
tic field in a superconducting microwave cavity [60|462| , 
the terminology and techniques of quantum optics will 
prove helpful in developing and explaining the model. 
The thermal reservoir is engineered by weakly coupling 
the harmonic oscillator with a succession of two-level sys- 
tems. By continually measuring their outgoing states, we 
are able to continuously monitor the evolution of the os- 
cillator. The effective equation of motion governing the 
dynamics of the continuously-monitored oscillator is a 
stochastic Schrodinger equation (Eq. [H] below), which is 
presented in this section. Its derivation, though a basic 
extension of the techniques utilized in Refs. p3L |63. [68|. 
is technical and is therefore reserved for Appendix [X] 

I have in mind a forced quantum harmonic oscillator of 
mass TO and frequency with position x and momentum 
p, coupled to a thermal reservoir at inverse temperature 
13, as depicted in Fig.[TJ The Hamiltonian of the harmonic 
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FIG. 1. Depiction of the energy levels E„{ft) for the instan- 
taneous eigenstates |n/j) of a forced quantum harmonic oscil- 
lator of frequency uj with forcing protocol ft. The oscillator 
interacts resonantly with a sequence of two-level systems with 
ground states |0) and excited states |1) in order to simulate a 
thermal reservoir at inverse temperature /?. 
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is parameterized by an externally controlled force, or pa- 
rameter, /, which is varied with time in order to do work 
on the oscillator. Here, and throughout the following 
h = 1. For an electromagnetic field, the forcing can be 
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induced by deterministically varying with time a collec- 
tion of macroscopic classical sources [tI] . As can be in- 
ferred from Eq. [1] a harmonic oscillator in the presence 
of an external force is equivalent to a harmonic oscillator 
that has been translated horizontally a distance /, and 
whose energy has been lowered by (l/2)ma;^/^. This ob- 
servation motivates the definition of lowering and raising 
operators parameterized by /, 
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(2) 



which when substituted into the Hamiltonian [Eq. [T] lead 
to the simplified expression 
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(3) 



For each /, the eigenstates \nf) of the number operator 
Nf = Oj^af, with eigenvalues n, are also eigenstates of 
the Hamiltonian H{f) with energies 



(4) 



From time i = to r, we drive the oscillator by varying 
the force / using the linear protocol = {ft}t=0' where 
ft = \/^j(rruJ)vt and v is the rate at which the force 
varies. For notational simplicity, I will denote the lower- 
ing and raising operators [Eq. [2] at time t when the force 



has the value /* as at = aj^ and 



= a^^. These opera- 
tors, at and aj, correspond to the lowering and raising op- 
erators for the eigenbasis \nf^) that at each t diagonalizes 
the Hamiltonian: H{ft)\nf^) ~ En{ft)\nft)- These states 
|n/j) - which are not solutions of the Schrodinger equa- 
tion - I will call the instantaneous eigenstates of H{ft), 
though_sometimes they are refereed to as the adiabatic 
basis ■'}{) due their employment in the derivation of the 
quantum adiabatic theorem. 

As we drive the oscillator, it will continually exchange 
energy with a thermal reservoir at inverse temperature 
/3. However, I will not model the thermal reservoir as 
very large number of degrees of freedom. Instead, I fol- 
low Refs. m, 113, lllj and engineer the thermal reservoir 
by weakly coupling the harmonic oscillator to a sequence 
of two-level systems, which I will call atoms since this 
procedure was originally proposed in the context of cav- 
ity quantum electrodynamics [60l - [6^ . A possible experi- 
mental realization in cavity quantum electrodynamics is 
illustrated in Fig. [5] The Hamiltonian of each atom is 



Ha ~ uja'^ (T, 



(5) 



where a — |0) (1| is the atomic lowering operator inducing 
transitions from the atom's excited state |1) to ground 
state |0). Prior to interacting with the oscillator, each 
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FIG. 2. Illustration of a possible experimental setup realiz- 
ing a forced harmonic oscillator weakly coupled to a thermal 
reservoir. A sequence of two level-atoms. A, pass through a 
cavity C, after which their final state is measured by detector 
D. 



atom is prepared in a definite known state: either its 
ground state |0) with probability 



1 



1 + e-'3" 

or its excited state |1) with probability 



1 + e-'3"' 



(6) 



(7) 



Then, one at a time, each atom interacts with the os- 
cillator for a short time St, after which we make an in- 
stantaneous projective measurement of the atom's state 
in order to determine whether it remained in its initial 
state or transitioned. By measuring the states of the 
atoms after they interact with the oscillator, we are able 
to track (or monitor) the oscillator's evolution. More- 
over, knowledge of whether the atom transitioned or not, 
allows us to monitor how much energy is transfered to 
the atom during its interaction with the oscillator, and 
consequently how much energy has been exchanged with 
the thermal reservoir. 

The coupling between each atom and the oscillator at 
time t is mediated by the displaced raising and lowering 
operator. 



Clf — Gf — iv jijj, 

through the interaction Hamiltonian 

i?/(/t) = A(al(T + atat), 



(8) 



(9) 



with weak-coupling, A ^ w. In the absence of 
driving (y = 0), this interaction Hamiltonian is the 
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Jaynes-Cummings Hamiltonian, known to quantitatively 
describe experimental observations in microwave cavi- 
ties [m [sOl- However, in the language of quantum op- 
tics, the lowering and raising operators of the harmonic 
oscillator are displaced by iu/uj due to the presence of 
the classical field induced by the motion of the classical 
sources. In addition, at and a\ are evaluated at t, in- 
ducing transitions between the instantaneous eigenstates 
of H{ft), which is vital for interpreting the sequence 
of atoms as a thermal reservoir constantly exchanging 
energy with the oscillator. A rigorous justification for 
the form of Hj in Eq. ^ is presented in Appendix |B] 
There a two-time scale perturbation analysis is applied 
to the Schrodinger equation for the oscillator coupled 
to an atom in the dipole approximation. The analysis 
indicates that the dynamics are well approximated by 
Hi{ft) [Eq. |9] on times XSt = 0(1) when the coupling is 
weak, A ^ a;. 

Having introduced the model, I now present the 
stochastic Schrodinger equation governing the evolution 
of the harmonic oscillator conditioned on the measure- 
ment outcomes js^, [56l - [59j . The stochastic Schrodinger 
equation (Eq. ll2l below) is an effective equation of motion 
for the evolution of the harmonic oscillator state vector 
{ipt) at time t conditioned on all the measurements made 
on the atoms prior to time t. In its derivation, I have as- 
sumed that (i) each atom interacts for the same amount 
of time St; (ii) during each St only one atom interacts with 
the oscillator; (iii) the interaction time is short, XSt ^ 1; 
(iv) the interaction is weak A <gC w; (v) the mean num- 
ber of excitations in the oscillator fit at each time t is 
small, XStfit <C 1; and (vi) the driving is not too fast, 
XSt{v/uj) ^ 1. While this last requirement restricts the 
rate of the driving i/, it does not necessitate the adia- 
batic approximation {v <^ uj): assumption (vi) can still 
be satisfied with v ^ lo. 

Consider the evolution of the harmonic oscillator dur- 
ing a small time interval dt large compared to St. During 
dt many atoms interact with the oscillator. Prior to in- 
teracting with the oscillator, each atom is prepared in 
a known state, either |0) or |1). As I demonstrate in 
Appendix \^ the probability during dt to measure an 
outgoing state different from an initial state ~ that is to 
observe a transition or jump - is small. In particular, the 
probability to observe a jump from |0) to |1) during dt 
is poi = 9^o{d\dt)tdt, where g = X^St is the jump rate 
and (•)( = (^/itl • \^t) denotes the quantum mechanical ex- 
pectation with respect to the conditioned oscillator state 
vector at t; the probability to jump from |1) to |0) during 
dt is pio — gri{dtd\)tdt] and the probability of observing 
no transitions is 1 — poi ~ Pw- Since the probability to 
observe a jump is of order dt^ the stochastic sequence of 
jumps is a Poisson process. Mathematically, I introduce 
two stochastic Poisson increments corresponding to an 
atom jumping up (the system jumping down), dA^+, and 
an atom jumping down (the system jumping up), dN^ . 
These Poisson increments are a sequence of random num- 



bers that are either or 1 - 

{dN+f=dN+ and {dN^f^dNt (10) 

- whose ensemble expectation value at t is 

E[dN+]^gro{d\dt)tdt, 
E[dN^] ^ gri{dta\)tdt. 



(11) 



The Poisson increments allow us to compactly write the 
stochastic Schrodinger equation for the stochastic differ- 
ential of the conditional state vector during the time in- 
terval \t,t + dt), d\'ipt) = IV't+dt) — IV't); as the Ito stochas- 
tic differential equation 

d\i>t) =dt(^ - iH{ft) ~ ^d\dt - ^ata\ 

, ff'^o ,_t- \ , ff^i/- -t\ Ai / \ 
+ —{alatit + —{atal)t j \ipt) 



dN: 




(12) 



Most of the time the stochastic increments dNl , j — 
+,— , in Eq. [1^ are zero and jV't) evolves deterministi- 
cally under the action of the term proportional to dt in 
Eq. [T21 Occasionally a jump is detected, signaling an 
abrupt change in the oscillator's state vector. When an 
atom jumps up dN^ equals one, and the oscillator's state 
vector experiences a discontinuous change mediated by 
at- Similarly, when dNt~ is one, the oscillator changes 
abruptly under the action of aj . 

The state vector |V't) characterizes the sub-ensemble 
of systems conditioned on a particular sequence of mea- 
surement outcomes. Because each measurement outcome 
is random, the conditioned sub-ensemble and its cor- 
responding state vector jV't) vary stochastically. Equa- 
tion [12] is the equation of motion describing this stochas- 
tic evolution of lipt) through Hilbert space [s^. The 
goal in formulating a quantum stochastic thermodynam- 
ics will be to associate with the stochastic evolution of 
the state vector and its corresponding ensemble notions 
of work, heat, and entropy. 

The state vector \tpt) describes the conditional (or se- 
lective) evolution of the oscillator. When the measure- 
ment outcomes are ignored (or averaged over), the uncon- 
ditional evolution is characterized by the density matrix 
obtained as the average over all measurements [59|, 



Pt^EUt){i^t 



(13) 



The equation of motion for the density matrix pt, calcu- 
lated by differentiating Eq. [13] with the aid of Eq. [12] is 
the Linblad master equation 59], 

dtpt = -i[H{ft),pt] + groV[at]pt + griV[a\]pt, (14) 
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where 



I?[c]/9 cpc^ — ■^^c^cp — -pc^ c. 



(15) 



For fixed / (j^ = 0), Eq.[T4]is a master equation describing 
the evolution of a harmonic oscillator in the presence of 
a thermal reservoir at the inverse temperature /3, whose 
value can be inferred from the ratio of Eqs. l6l and [71 [59| . 



(16) 



In particular, if we allow the system to evolve freely with 
/ fixed, in the long-time limit the system will relax to 
equilibrium with a density matrix given by the Boltz- 
mann density matrix 



en 
Pf 



'l3H(f) 



(17) 



where Zf = Tr[exp(— /3_ff(/))] is the partition function. 
Thus, the average evolution of the harmonic oscillator 
induced by its interaction with a sequence of two-level 
atoms is equivalent to the evolution caused by a thermal 
reservoir. In general, when v ^ Q, Eq. [14] describes the 
nonequilibrium evolution of the oscillator. Specifically, 
the solution to the master equation [Eq. [TJ] at t, pt, may 
differ from the instantaneous equilibrium density matrix 

en 
Pf.- 



III. QUANTUM STOCHASTIC 
THERMODYNAMICS 



In this section, I introduce a quantum stochastic ther- 
modynamics for the forced harmonic oscilator outlined 
in Sec. [Ill Stochastic thermodynamics, like its classical 
counterpart, will be formulated along thermodynamic 
processes, which are introduced in Sec. IIII Al Then in 
Sec. [IITbI the change in energy along a thermodynamic 
process is divided into two contributions: the work and 
the heat. Entropy and entropy production are identified 
in Sec. [mCl 



A. Thermodynamic Processes 

In classical macroscopic thermodynamics, work and 
heat characterize the exchange of energy between a sys- 
tem and its surroundings during a thermodynamic pro- 
cess. Such a process is specified by a sequence of macro- 
scopic actions (or events) executed to manipulate the 
evolution of the system. Typical actions taken during 
a thermodynamic process include preparing the system 
initially in thermodynamic equilibrium at fixed temper- 
ature; or modifying the volume or the pressure using 
a specified protocol. Similarly, in order to construct a 
quantum stochastic thermodynamics, I will first specify 



the sequence of actions that define a thermodynamic pro- 
cess. Two processes, introduced below, will prove useful: 
the forward and reverse processes, which are related by 
time reversal. 

A process will be composed of four parts: an initial 
preparation of the system; an initial observable that is 
measured at the beginning of the process; a protocol for 
varying the external force with time; and a final observ- 
able measured upon completion of the process. The mea- 
surements of observables at the beginning and end of the 
process are macroscopic actions unique to quantum pro- 
cesses. They are necessary in order to guarantee that 
the forward and reverse processes begin in a pure state, 
and are required in the proof of the detailed fluctuation 
theorem in Sec. IIVI below. 

We begin the forward process by preparing the initial 
ensemble of quantum systems. We form this ensemble 
by collecting a large number of quantum systems each 
in a known eigenstate |a) of an observable A with non- 
degenerate eigenvalue a, such that the fraction of sys- 
tems in eigenstate \a) is proportional to the probability 
Pa- The statistical properties of this ensemble are char- 
acterized by the density matrix 



Po 



= PA = ^Pa\a){a\. 



(18) 



We then randomly select a quantum system from the 
initial ensemble pA, which will be in state \a) with prob- 
ability Pa- To confirm this state, we make a projective 
measurement of A. Then from time t ~ to t, the force 
is varied according to the protocol while continuously 
monitoring the oscillator through the sequence of inter- 
acting atoms. During this time the oscillator state vector 
evolves stochastically according to Eq. [T^l In any real- 
ization, we will observe that the atoms jump at a specific 
sequence of times. We record each jump time tk and the 
type of jump nik = — , where -I- denotes an atom jump- 
ing up and — denotes an atom jumping down. At time 
T, we fix the external force at fr — \J 2 j (miij)vT , cease 
the flow of atoms, and make a projective measurement 
of a new observable B. The outcome is one of its non- 
degenerate eigenvalues 6, corresponding to the eigenstate 
1 6), with probability P^ ~ . Consequently, at 

time t = T the state vector for any realization will col- 
lapse into one of the eigenstates of B, and the density 
matrix [Eq. [13] will be diagonal in the eigenbasis of B: 



PB=J2Pb\b){b\- 



(19) 



Repeating this series of actions generates an ensemble of 
realizations. 

During each realization of the process, we make a mea- 
surement of A with outcome a at the beginning of the 
process, a measurement of B with outcome b at the 
end of the process, and observe a series of jumps ruk at 
times tk- I collect this string of measurement outcomes 
into a vector, which I call the measurement trajectory. 
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7 = {a;mi,ti;m2,t2; ■ ■ ■ ;^}- Furthermore, at each time 
t the system is characterized by the state vector condi- 
tioned on all measurements up to t, 1^4(7)), which is the 
solution of Eq. [T^ The quantum trajectory is the se- 
quence of state vectors traced out by the system through 
Hilbert system during the time interval from t = to t, 
denoted by ipil) = {\4't{j))}t=o '^^ simply -0 when 7 is 
clear from the context. 

The reverse process is defined as the time-reversed for- 
ward process in which each action of the forward pro- 
cess is carried out in the reverse order. First, recall that 
time-reversal in quantum mechanics is implemented by 
the time- reversal operator Q, which is an antilinear - 

ec = c*e, (20) 

for any complex number c, where * denotes complex con- 
jugation - involution, 

= /, (21) 

where I is the identity operator [tt} . 

In the reverse process, we prepare the initial ensem- 
ble by collecting a number of quantum systems, each 
in an eigenstate of the time-reveresed observable of B, 
B = QBQ~^. The ensemble is constructed so that each 
eigenstate \b) — Q\b) oi B with eigenvalue b occurs with 
probability Pb ~ which in general differs from the prob- 
ability to measure b at the end of the forward process, 
Pb- Hence, the initial ensemble is characterized by the 
density matrix 

Po = PB=^n|&)(&|. (22) 

6 

We then randomly select a quantum system from this 
ensemble and make a projective measurement of B. Af- 
ter, we drive the system from time t — to t using the 
time-reversed protocol T = {ft}t=o^ where ft = fr-t — 
^Ij (jruxi)v(T — t). During this time interval, the oscil- 
lator interacts with a sequence of atoms; we record the 
times at which the atoms jump and the types of jumps 
mfc, as in the forward process. Finally, we measure the 
time-reversed operator A = QAQ~^ obtaining eigenvalue 
a corresponding to eigenvector \a) = 8|a) with probabil- 
ity Pa = |(a|V'r-)P- At the time of measurement the 
state vector collapse into an eigenvector of A and the 
density matrix becomes diagonal in the eigenbasis of A: 

a 

Every measurement trajectory of the forward pro- 
cess during which M jumps are observed, 7 — 
{a; mi, ti] m2, ^2; • • • ; 'mM,tM] is paired with a conju- 
gate reverse measurement trajectory of the reverse pro- 
cess 7 = {b\Thi,T-tM; ■ ■ ■ ;mM-i,T-t2;rhM,T-ti;d}, 
where iriM-k is + (— ) when m^+i is — (-f), for k — 
0, . . . M — 1: when an atom jumps up along a forward 



trajectory, the atom will jump down in the conjugate 
reverse trajectory. The quantum trajectory traced out 
through Hilbert space by the oscillator's state vector cor- 
responding to the measurement trajectory 7 will be de- 
noted ■(/'(t') — {\^t{l))}t=oi where |V't(7)) is the solution 
of Eq. [12] during the reverse process. Observe that in 
general the sequence of state vectors traced out by the 
system is not the time-reversal of the quantum trajectory 
observed in the forward process, |V't(7)) 7^ 0|V'r-t(7))- 

B. First Law of Stochastic Thermodynamics 

In this section, I formulate a energy balance equation 
(Eq. 133] below) relating the heat and work to the change 
in energy along individual quantum trajectories. I pro- 
ceed by first identifying the heat as the energy transferred 
to the thermal reservoir. The energy is then defined as 
the quantum mechanical ensemble average of the Hamil- 
tonian. Finally, the change in the energy not accounted 
for by heat is identified as the work. 

Heat is the energy exchanged with a thermal reser- 
voir. Since the thermal reservoir in the present model 
is composed of a sequence of two-level atoms, I identify 
the heat absorbed by the oscillator as the energy released 
from the atoms. Energy is only exchanged when an atom 
transitions or jumps between its energy eigenstates. In 
particular, each time we observe a jump from |0) to |1), 
the atom has absorbed uj energy from the oscillator; sim- 
ilarly, when the atom has jumped from |1) to |0), it has 
released w energy. Therefore, for a given quantum tra- 
jectory '0(7)1 the increment in the heat absorbed by the 
oscillator during the small time interval from t io t + dt 
may be written as 

dQt[ilj] = i^{dNt ~ dN+). (24) 

The heat absorbed by the oscillator during [0,t) is the 
stochastic integral of Eq. [24l 

Qt[^] = / dQs[^] (25) 

JQ 

= f Lo{dN--dN+). (26) 
Jo 

Next, I define the (internal) energy at time t for the 
entire ensemble characterized by the state vector |'0t(7)) 
conditioned on all the measurements up to The energy 
at t along the quantum trajectory 0(7) is the ensemble 
average of the Hamiltonian, 

= {ijt\H{ft)\iPt). (27) 

The energy varies with time due to the time dependence 
of the Hamiltonian as well as the stochastic evolution of 
the state vector. In particular, the change in the energy 
up to t is 

AC/t[V] = / dUsi^j] (28) 
Jo 
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where, dUt[tp] = Ut+dtli^] — Ut[tp], is the stochastic differ- 
ential of Ut during the small interval from t to t + dt: 



dUtiij] ^dtft{HdfH{ftMt 



dtgro{ipt\ 
dt gri {ipt I 



^{H{f,),alat}-H{ft){alat)t 
l-{H{ft),atdl}~H(ft){atdl)t 




atH{ft)dl)t 

{dt4)t 
d\H{ft)-at)t 




(29) 



which can be deduced from the definition of Ut [Eq. [28] 
with the aid of Eq. [T^l Here, for two operators O and 
O', {O, O'} = OO' + O'O is the anti-commutator. For a 
process where the initial observable is the initial Hamil- 
tonian, A = H{fo), and at the end of the process the 
final Hamiltonian is measured B = H^fr), this defini- 
tion of the change in internal energy in Eq.[5S] agrees with 
that commonly encountered in the derivation of quantum 
work relations [l^. However, its statistics are diferrent 
due to the interaction with the reservoir. 

Finally, the work during the time interval from t to t + 
dt is identified as the change in the energy not accounted 
for by heat 



dWtii^] = dUt[^] - dQt[^], 



(30) 



which may be formulated, using the definitions of 
dUt [Eq.Hg and dQt [Eq.lM], as 



dWtli^] ^dtft{MdfH{ft)\iPt 



dtgrQ{ipt\ 
dt gri {ipt I 



^{H{ft),dlat}-H{ft){dldt)t 
l-{H{ft),at4}-Hift){dtal)t 



■dN; 



-dN+ 



{atHUt)d\) 
{dta\)t 

{4H{ft)dt) 
{aldt)t 



{H[ft))t 
mh))t 



(31) 



The work along a quantum trajectory during the time 
interval [0,i) is the stochastic integral of dWt, 



(32) 
(33) 



Jo 



Remarkably, the average of Eq. [32] over all realizations, 
after substituting in Egs. [TT]and[5Tl has the simple form 

E[Wt]^ f dthE[{i^t\dfH{U)\i,t)l (34) 
Jo 



reminiscent of the definition of work for classical sys- 
tems. Equation [M] resolves an ambiguity in a common 
definition of the average work along nonequilibrium pro- 
cesses [ts'I. One method for identifying the work is to 
differentiate the average energy ut — TT[H{ft)pt] with 
time, 

dtut = Tr[dtH{ft)pt] + TT[H{ft)dtPt], (35) 

and then to identify the average work as wt — 
Tr[dtH{ft)pt] and the average heat as qt ^ Ti[H{ft)dtPt]- 
However, when the trace is evaluated using a time- 
dependent basis, Eq. [35] is no longer correct and the 
seperation of dtUt into two parts is not unique [47| . Equa- 
tion [34] avoids this ambiguity, as it is evaluated using 
state vectors instead of density matrices. 

To illustrate the relationship between work, heat, 
and energy, I have numerically integrated the stochas- 
tic Schrodinger equation [Eq. [12] for a thermodynamic 
process where the initial energy A = H{fo) and fi- 
nal energy B = H{fr) are measured [HI. Plotted in 
Fig. [3] is the change in energy Al/j [Eq. [55], heat Qt 
[Eq. [25], and work Wt [Eq. [32] as a function of time 
from i = to r = 80s along a representative quan- 
tum trajectory beginning in energy eigenstate \'2.fg) and 
found at r in state The work to lower the energy 

zero, Wt = —{l/2)muj'^j^dsdsfs ~ —ujv'^t^ , has been 
neglected in Fig. [3] in order to highlight quantum effects. 
At two times, gt « 2.6, and 4.9, there is a discontinuous 
change in the average energy in Fig. [3l^a) when an atom 
jumps and exchanges a quantum of heat with the reser- 
voir, as depicted in Fig. [3l^6). In between the jumps, the 
average energy varies rapidly taking on non-integer val- 
ues, a signature that the oscillator's state is in a coherent 
superposition of instantaneous energy eigenstates. This 
additional variation of the energy is accounted for by the 
work in Fig. [3l^c). Furthermore, when there are discon- 
tinuous jumps the change in the energy is not completely 
due to the flow of heat. As can be seen in Fig. [HJc), the 
work Wt changes discontinuously as well. These jumps 
are a consequence of the oscillator being in a superposi- 
tion of energy eigenstates. However, the additional work 
accrued during these jumps is not work due to the varia- 
tion of the external parameter ft , as is typical of classical 
systems. Its origin is the sudden change in the state vec- 
tor triggered by the measurement of the outgoing atom, 
though how this energy is transfered to the oscillator 
remains unclear. Nevertheless, this definition of work 
provides a consistent framework to describe the flow of 
energy between the oscillator and its environment. 



C. Second Law of Stochastic Thermodynamics 

The second component of a quantum stochastic ther- 
modynamics is a definition of entropy and entropy pro- 
duction along individual quantum trajectories. 

I first define the change in entropy of the thermal reser- 
voir. The role of the thermal reservoir is played by the 
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as heat, 




FIG. 3. Plots of (a) change in energy A?7t, (6) heat Qt, and 
(c) work Wt as a function of time t from f = to r = 80s along 
a represen tative quantum trajectory with force protocol ft = 
(rnijj)vt initially in energy eigenstate I2/13), whose final 
state is measured to be 1 1/^ ) . The additional work required to 
lower the energy zero, wt = —lov^t^, has been neglected. The 
parameter values are uj = v = IQ Hz, A = 0.5 Hz, St = 0.1s, 
g — 0.025 Hz, ri/ro — 0.75, and m = 1kg. 



sequence of atoms; they act as a very large depository 
for energy with inverse temperature /3. Specifically, be- 
fore each thermodynamic process we prepare a large col- 
lection of atoms each either in their ground or excited 
states with relative probability given by the Boltzmann 
weight [Eq.[T6]. Over the course of the process, each atom 
interacts with the harmonic oscillator, but the fraction 
of atoms that jump is small. As a result, the relative 
fraction of atoms in the ground and excited states will 
not deviate appreciably from the Boltzmann distribution. 
Consequently, the change in entropy of the reservoir at 
time t along a particular quantum trajectory "0(7) is pro- 
portional to the energy absorbed by the thermal reservoir 



dsim^-pdQtm. 



(36) 



Integrating gives the total change in reservoir entropy 
along the process 



dsl[^]=-pQrm- 



(37) 



In classical stochastic thermodynamics, the entropy of 
the system is associated with the Shannon entropy of the 
system's phase space density. The quantum version of 
the Shannon entropy is the von Neumann entropy, which 
at time t reads 



Sipt 



-Tr[pt In pt] 



(38) 



From the definition of the density matrix in Eq. 1131 we 
may rewrite the von Neumann entropy as a classical sta- 
tistical average over quantum trajectories V'(7) 



S{pt)^-E[{^Jt\lr^Pt\iJt)l 



(39) 



where here pt is the density matrix for an ensemble of re- 
alizations [Eq. I13| and evolves deterministically accord- 
ing to the master equation in Eq. 1141 The form of S in 
Eq. [31] sug gest s defining the trajectory-depedent system 
entropy as [73| 



St[ip] ^ -{ipt\ InptlV't 



(40) 



Following Seifert 0, one may attempt to develop a 
stochastic differential equation for st (as in Eq. for 
dUt); however, there is no obvious compact expression, 
since pt in general does not commute with its time deriva- 
tive. Over the course of a trajectory that begins in state 
\a) and ends in state \b) the change in the system's en- 
tropy is 



As[i)] = -(■0t| In Pr I ^/^r) 

= -lnPfc + lnP„. 



(■00 1 In Po|0o) 



(41) 



where in the second line I used that the density matrix is 
diagonal at the beginning [Eq. (TB] and the end [Eq. [TO] 
of the process. 

Adding Eqs. [37land[4Tl we find that the total entropy 
production during a thermodynamic process is 



Astot[ip] = As[ij] + Asrlip] 

= - In Pb + \n Pa ~/3Qr[iP], 



(42) 



consistent with the definition proposed by Monnai [80| 
In Sec. lIVl I show that the average total entropy produc- 
tion AStot — E[Astot] is a measure of the irreversibility 
of the thermodynamic process and is non-negative in ac- 
cordance with the second law of thermodynamics. Fur- 
thermore, the average total entropy production agrees 
with the definition introduced and analyzed previously 
by Breuer [8l|. 



9 



IV. DETAILED FLUCTUATION THEOREM 

In this section, I address the detailed fluctuation the- 
orem [2, [ll [13, [13, ii-ii] in the context of quantum 
trajectories. The detailed fluctuation theorem relates 
the probabilities to observe particular microscopic tra- 
jectories along two thermodynamic processes related by 
time-reversal. It has been derived for a wide class of dy- 
namics; in each case, the detailed fluctuation theorem 
identifies the source of time-reversal symmetry breaking 
during a thermodynamic process with the total entropy 
production [s^, Hs, H^l- The generality of this obser- 
vation suggests that we may use the detailed fluctuation 
theorem as a tool for ascertaining the entropy production 
during a nonequilibrium thermodynamic process in situ- 
ations where its definition may not be obvious. Through 
verifying the detailed fluctuation theorem for quantum 
trajectories, we will see that the source of time-reversal 
symmetry breaking is the same as the total entropy pro- 
duction Astot [Eq-112]; thereby, providing evidence that 
the definition of entropy production in Sec. IIII CI is con- 
sistent with its role as a measure of irreversibility. 

The detailed fluctuation theorem for the quantum tra- 
jectories of our forced harmonic oscillator relates the 
probability to observe a quantum trajectory '0(7) in the 
forward process P[4'] to the probability to observe the 
conjugate trajectory "0(7) in the reverse process P[V'] as 

lnJf^ = As,otm, (43) 

derived below. In quantum systems, the detailed fluctu- 
ation theorem has been demonstrated previously in nu- 
merous situations g, [H [H [M IHil IH, Ho]. The 
essential ingredients are the time-reversal symmetry of 
the Hamiltonian, and that projective measurements are 
made at the beginning and at the end of the pro- 
cess. Moreover, any collection of measurements may 
be performed during the process without invalidating 
Eq. HSlliil. In this respect Eq. |33]is not novel; since, the 
Hamiltonian in Eq.[T]is time- reversal invariant; a thermo- 
dynamic process, by definition, begins with the measure- 
ment of the observable A and is terminated by a mea- 
surement of the observable B; and the conditional evolu- 
tion of the state vector \tpt) embodied by the stochastic 
Schrodinger equation is merely a sequence of weak mea- 
surements. Nevertheless, I will sketch the derivation of 
Eq. 1431 since alternative approaches sharpen our under- 
standing. 

Before deriving Eq. 1431 let me comment on one of its 
consequences. In particular, notice that Eq. immedi- 
ately implies that AStot — E[Astot] equals the relative 
entropy Dif\\g) = J dxf{x) \n[f{x)/g{x)] of P and P, 

AStot = DiP\\P) > 0, (44) 

which is always non-negative {D > 0) [9l|. The relative 
entropy D{P\\P) is a measure of the distinguishability of 



P and P. Therefore, AStot is a measure of how distin- 
guishable a forward process is from its time reverse '93] ; 
it measures the thermodynamic irreversibility of a pro- 
cess [90*1. In particular, the equality in Eq. |31]is reached 
only for reversible processes where the forward process is 
indistinguishable from the reverse process, P[i/'] — P[ijj]- 
Note, unlike the entropy in macroscopic thermodynam- 
ics, which enters into the second law of thermodynam- 
ics and also encodes information about the properties 
of equilibrium states such as specific heats, the entropy 
defined here only reflects the irreversibility of a nonequi- 
librium thermodynamic process; except when the ther- 
modynamic process begins and ends in equilibrium. 

To verify the detailed fluctuation theorem [Eq. 031 , we 
calculate P[ip] and P[4>]- To determine P and P, let 
me first note that within the quantum trajectory for- 
malism the stochastic evolution of the state vector is 
the result of a sequence of weak measurements. Each 
measurement outcome can be represented by measure- 
ment operator [9^, whose effect on any pure state 
jf/;) leads to an unnormalized state vector \'^) = 
that encodes the probability for that outcome in its norm 
p(ri) = (\E'|^) = {ijjlQ'f fl\tp) . Similarly, for a sequence of 
measurements, the norm of the unnormalized state vec- 
tor obtained by applying a series of measurement oper- 
ators equals the probability to observe that sequence of 
outcomes. Therefore, the probability to observe any so- 
lution of the stochastic Schrodinger equation, P[ip], can 
be found by determining the norm of the unnormalized 
state vector resulting from the action of a series of mea- 
surement operators describing the effects of observing a 
series of jumps or no jumps in the sequence of atoms. 
These measurement operators can be obtained from the 
structure of Eq. [12] jsl, 113, [13, Ell . Below, I simply re- 
port the results in order to keep the present discussion 
concise. 

In the forward process, the evolution of the unnormal- 
ized state vector in between jumps is realized by an ef- 
fective time-evolution operator Ueff{t,s) - the solution 
of 

dtU,ffit,s) = -zHe//(/t)t^e//(i,s), (45) 

with non-hermitian effective Hamiltonian 

H^ffif) = Hif) i^a\af - z^afa\ (46) 

with initial condition U (s, s) = /. This deterministic 
evolution is punctuated by discontinuous changes when 
a measured atom jumps, induced by the jump operators 

j- (t) = Vdtgnal ^^^^ 
j+{t) = \/dt groat, 

which satisfy 

Ut)=jl{t)eP-/'. (48) 

Notably, Eq. |351 originally derived by Crooks in a 
more general setting assuming the thermal reservoir is 
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in equilibrium, may be viewed as an operator extension QH{f)Q~^ , implies the relationship between effective 
of detailed balance to quantum trajectories, relating an time-evolution operators, 
atom's transition between energy eigenstates to the en- 



ergy absorbed by the thermal reservoir as heat. With 
this notation, the unnormalized state vector at time r 
conditioned on 7 = {a; mi, ti; m2, ^2; • • • ]TnM , tM] 6} is 



\^r{l)) =\b){b\Ueff{T,tN)j^AtM) 

^\b){b\C[^)\a), 



and the probability to observe 7 given that the initial 
state is \a) is the norm of 



F[^|a] = (VI., (7) I (7)). 



(50) 



In the reverse process, the force protocol is reversed, 

implying that ft = —fr-t oc —v and that each appear- 
ance of v in the forward process must be replaced by — 
in the reverse process. Consequently, the reverse effective 
Hamiltonian at time t that generates the deterministic 
(but non-hermitian) evolution Uef f{t, s), is 



.rig. 



(51) 



where the reverse displaced raising and lowering opera- 
tors are obtained from Eq. [8]by the substitution v ■(-> —v 
as 



at 



+ iv jijj. 



(52) 



The discontinuous jumps are induced by reverse jump 
operators 



j- (t) = y^dtgnal 
j+{t) = yjdtgrodt, 

related to the forward jump operators by 

j_(t) = ej|(r-oe-V-/^ 



(53) 



(54) 



C/e//(t,s)=eC/J^/7 



(58) 



■JmAtl)Ueff{tlMa) 

(49) 



Combined with Eq. [M] this leads to a time-reversal sym- 
metry between C [Eq. |49] and C [Eq. [56] 



The probability to observe the conjugate quantum tra- 
jectory ^"(7) conditioned on initiating the trajectory in 
state \b) is given by the norm of the unnormalized state 
vector 



£(7) = e£t(^)e-ie- 



(59) 



which when substituted into the definitions of P [Eq. [ST 
and P [Eq. [57] , gives 



ln:^i^ 



P[^\b] 



(60) 



Recalling that in a thermodynamic process the proba- 
bility for a quantum trajectory to start in state \a) is 
Pa [Eq. [H] and the probability in the reverse process to 
begin in state |6) is Pb [Eq. [22], we find that the log of 
the ratio of P[-!/;] ^ P[i^\a]Pa and P[-0] ^ P[^\b]Pb is 



(61) 



Equation [61] is valid for any Pa and Pb. However, as in 
classical stochastic thermodynamics, there are two note- 
worthy special cases |2j]. First, when the initial density 
matrix of the reverse process is the final density matrix 
of the forward process, po = Pt = Pb, we recover Eq. [43] 
equating the total entropy production to the irreversibil- 
ity of the process. A second special case of Eq. [51] con- 
nects the work dissipated to irreversibility. Consider a 
forward process where the initial observable is the initial 
Hamiltonian, A = H{fo), the final observable is the final 
Hamiltonian B = H(fT-), and both the forward and re- 
verse processes are started in equilbirium: pQ = p^^ and 
Po — pj' ■ In which case, Eq. [HI] reduces to 



lnW=/3(iy.[^] 

P[v\ 



AF), 



(62) 



\a){a\Cmb), 



(55) 
(56) 



as 



P[m = (^'.(7)l*r(7))- 



(57) 



where /S.F is the equilibrium free energy difference be- 
tween the equilibrium ensembles with external parameter 
values /o and Thus, the work dissipated - the work 
done in excess of the free energy difference - has a clear 
physical interpretation as a measure of the irreversibility 
iAr,0^£)^ thermodynamic process where the system is driven 
between two equilibrium states. 

To conclude this section, I comment on the work fluc- 
tuation relations. In particular. Eg . 162 1 immediatelv leads 
to a work fluctuation relation relating the probability 
to observe W work in the forward process p{W) to the 
probability to observe —W work in the reverse process 
p{-W) HI 



To complete the derivation of Eq. [33] we note that 
the time reversal invariance of the Hamiltonian, H{f) = 



p{W) 
p{-W) 



^I3{W-AF) 



(63) 
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and its integral version, the nonequilibrium work fluctu- 
ation theorem [21], 

(e-'3(w^-^^)) = 1, (64) 

where here the angled brackets (•) denote an ensemble 
average over work values. Similar considerations apply to 
a detailed and integral fluctuation relation for the total 
entropy production 

V. DISCUSSION AND PERSPECTIVES 

In this article, I have formulated a quantum stochastic 
thermodynamics for quantum trajectories. As a concrete 
example, I investigated a forced quantum harmonic oscil- 
lator. The oscillator was coupled to a thermal reservoir 
composed of sequence of two-level atoms, which allowed 
us to monitor the energy transfer to the atoms and con- 
sequently the thermal reservoir leading to a physically 
motivated definition of heat. In essence, we are repeat- 
edly measuring the environment and not the system di- 
rectly. However, since the environment is broken up into 
individual quantum systems that interact with the oscil- 
lator one at a time, we are able to avoid any practical and 
conceptual difficulties with measuring an infinite thermal 
reservoir. Following the definition of heat, I introduced 
path-dependent definitions of work and entropy for indi- 
vidual quantum trajectories in Sec. IIIII Their connection 
with irreversibility was made in Sec. IIVI through the de- 
velopment of a detailed fluctuation theorem for quantum 
trajectories. 

The present discussion focused on a particular model. 
Nevertheless, the formulas presented in Sec. IIIII for work, 
heat, and entropy could be applied to other systems 
within the quantum trajectory formalism. The cen- 
tral issue is whether the measurement scheme can be 
given a consistent thermodynamic interpretation, as in 
the present model where the atoms played the dual 
role of quantum probe and thermal reservoir. This 
opens the interesting question of what other monitor- 
ing schemes (so-called unravelings of the master equa- 
tion), have consistent stochastic-thermodynamic inter- 
pretations m, [13, m HI . 

Within the quantum trajectory formalism, quantum 
trajectories are the stochastic evolution of the system's 
state vector through Hilbert space. A seemingly alterna- 
tive deflnition of a trajectory is offered by the consistent 
(or decoherent) histories 95^ interpretation of quantum 
mechanics. However, Brun has demonstrated that quan- 
tum trajectories can be consistentl y re presented within 
the consistent histories framework [96i|. A distinct per- 
spective is obtained by viewing these quantum trajecto- 
ries as a particular unitary dilation of the damped har- 
monic oscillator (97j . A unitary dilation of a quantum 
Markov semigroup is a representation of the irreversible 
dynamics as a unitary evolution on a larger Hilbert space, 
such that when projected onto the system's Hilbert space 



we recover the original irreversible evolution. Here, the 
irreversible dynamics of the harmonic oscillator is dilated 
onto the unitary dynamics of the oscillator plus the se- 
quence of atoms (see for example Ref. [695). 

A principle motivation for the present study was to de- 
velop trajectory-dependent deflnitions of work, heat, and 
entropy that could be addressed experimentally. The 
quantum harmonic oscillator can be experimentally re- 
alized as a single mode of an electromagnetic field in 
a microwave cavity (60l-l62j. and the thermal bath can 
be engineered by passing a sequence of two-level atoms 
through the cavity one at a time. In particular, the 
work fluctuation relations can be verified experimen- 
tally, since it is feasible to construct and measure with 
a quantum non-demolition measurement the individual 
energy eigenstatcs of the electromagnetic field ^60]. The 
greatest difficulty in verifying the predictions of quan- 
tum stochastic thermodynamics is the efficient detection 
of the atoms once they have interacted with the field. 
Quantum stochastic thermodynamics requires a near per- 
fect detection efficiency; whereas, modern experimental 
setups only reach forty percent [6^. Nevertheless, quan- 
tum stochastic thermodynamics could be investigated in 
other quantum systems where the quantum trajec tory 
formalism has been applied, such as quantum dots |98| . 
nanomechanical resonators f9^, or perhaps in quantum 
cyclotron s - w here thermal quantum jumps have been 
observed [lO0j |. 

Future research directions are manifold. Particularly 
interesting is a thermodynamic analysis of effects with 
a purely quantum origin; such as, the work required 
to generate entanglement, or the entropy produced in a 
squeezed thermal bath. The stochastic Schrodinger equa- 
tion in Eq. [12] for the linearly- forced harmonic oscillator 
may be of interest in its own right. For example, in the 
thermodynamic adiabatic limit, Eq.[T2] could be used to 
analyze geometric phases in open quantum systems. 
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Appendix A: Derivation of the Stochastic 
Schrodinger Equation 

In this appendix, I derive the stochastic Schrodinger 
equation [Eq. [l^] for the time evolution of the one- 
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dimensional forced harmonic oscillator conditioned on 
continuously monitoring its interaction with a sequence 
of two-level a toms , by adapting the method used in 
Refs. [H, [sl, \lQl\ . The stochastic Schrodinger equa- 
tion is an effective equation of motion that is valid for 
times long compared to the oscillator-atom interaction 
time St; given that only one atom interacts with the har- 
monic oscillator at a time; each atom interacts for the 
same amount of time, 5t; the interaction time is short, 
XSt <^ 1; the coupling is weak, A ^ w; and the mean 

number fit = {a\at)t is small XStfit <C 1. 

The derivation of the stochastic Schrodinger equation 
proceeds by evolving the oscillator-atom state vector, \xt) 
for a short time from <o to + 5t using the Schrodinger 
equation. We then calculate the probabilities to measure 
the atom in its ground and excited states in order to 
demonstrate that the atom's transitions are described by 
a Poisson process. 

The oscillator-atom state vector at time i, \xt)^ is the 
solution to the Schrodinger equation 



^^'H{ft)\xt), 



(Al) 



dt\xt) 

where the Hamiltonian 

= uj ^a^fl/ + — ^ — — mw^/^+w(T^(T+A(ay(T-f a/cr^) 

(A2) 

is a sum oi H{J) [Eq.[3], Ha [Eq.[5], and Hi{J) [Eq.[9]. 
Since each atom is prepared independently from the os- 
cillator, the initial condition at to is a factorized state, 
Ixto) = \'4'io)\4>ta) ^ with harmonic oscillator state vector 
l^/ito) and atom state vector \4>to)- 

The analysis is facilitated by first switching to the 
adiabatic-interaction picture jl02l |. To this end, let me 
introduce the translation operator 

A{f) = e-^^f, (A3) 

whose action on the position operator x and lowering 
operator is given by 

A\f)xA{f)^x + f, 

A\f)afA{f) = -a. 

Using A, I introduce the adiabatic-picture state vector 

\xt}A=A^{ft)\xt}, (A5) 

whose equation of motion is found by substituting Eq. lAll 
into the time derivative of \xt)A [Eg. IA5| and exploiting 
the properties of A in Eq. IA4[ 

dt\xt)A^~tHA{ft)\xt)A, (A6) 

where the adiabatic Hamiltonian is 

HAift) =A\m{ft)A{ft) + tA^{ft)A{ft) (A7) 



(A4) 



=LU [ a'a + i j - ^m.Lo'^ff + Loa^a (A8) 

^ pft + A(aV + aa^). 

Next, we shift to the interaction picture using the oper- 
ator K{t,to), defined as the solution to the differential 



equation 

dtK{t,to) = 



2 _f 2 
t 



u) ( a^a-\ — I muj^ f, 

\ 2/2 

+ WCrV -pf, 



Kit, to). 

(A9) 



with initial condition Kjtn , tp) = /. The solution may be 
obtained analytically (l03j | 



K{t,to) =exp 



-iu! a'a + (tV + 7: ] (t — to) + - 



dxfl 



to 



where, recalling that ft = ■\/2/ {rrujj)vt, 
dxf^e"^^='-*°'> 



= I— (e 
CO V 



2 .'to 



1 



(AlO) 

(All) 
(A12) 



Pt.ta — ^ 



dx dy fyfxsm[uj{y ~ x)] (A13) 

to "'to 

— \t-to--sm[ij{t-to)]\, (A14) 



and 



D{r]) — exp{T]a^ — 77*0) 



(A15) 



is the displacement operator, whose effect on the lowering 
operator a is 

D^{ri)aD{ri) = a + rj. (A16) 

The adiabatic-interaction state vector is defined as 

\Xt)Ai = K\tM)\xt)A. (A17) 

The equation of motion for \xt)Ai is obtained by differ- 
entiating with time Eq. IA171 and then substituting in 
Eqs. EH ESl EH Ell followed by EMI to conclude 
that 



dtWtlAi = -'iHAi{ft)\xt}Ai, 
where the adiabatic-interaction Hamiltonian is 



(A18) 



HAiift) = K^t, to)HA{ft)K(t, to) + ikHt, to)K{t, to) 

(A19) 

= A(atcr-f acr^). (A20) 



The next step is to solve Eq. IA18I perturbatively for 
short times. To clearly identify the approximations in- 
volved, I scale time s = At; scale the adiabatic-interaction 
Hamiltonian Hai = Hai/X; and scale the force, ft = 
{muj)gut, by introducing the dimensionless function 
g^t — i^t- In terms of these scaled quantities, the formal 
solution of Eq. IA18I - suppressing the subscripts AI in 
order to simplify the notation - is 
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fc=l ^0 



dskH{gps^)---H{gps^)\Xs 



(A21) 



where v = v/uj. Setting s = sq + <5s and expanding Eg. IA2ll to second order in 5s = XSt <C 1, we find 



Ss^ 

~ \Xso) - iSs{a''a + aa'')\xso) ~ —ia''''^ + aa'')'^\xso) 



(A22) 
(A23) 



where in the second hne I substituted in the definition of 
Hai in Eq. IA201 The validity of the asymptotic expan- 
sion in Eq. lA22l for Ea. lA23p requires that each successive 
term is smaller than the previous. In particular, we must 
demand that 

6s- < 1. (A24) 

Since 6s <C 1, i^/uj cannot be too large, restricting the 
rate at which the force varies. 

Now, let us consider the case when the atom is initially 
in its ground state at sq. Substituting |0S(,) = |0) into 
Eq. [M3l leads to 

\Xso+Ss) -\Tpso)\0) - i6sa\^so)\^) - — a^a|V'so)|0). 

(A25) 

We then measure the state of the atom. After the mea- 
surement, the unnormalized state vector of the harmonic 
oscillator conditioned on the outcome |0) is obtained by 
projecting Eq. I A2 5 1 onto |0), 

IM^o + 6s)) ^ - ^-^a^d\^so) = ^oHso)- (A26) 

Likewise, when the atom is found to be in the excited 
state, the unnormalized state vector of the harmonic os- 
cillator conditioned on the outcome 1 1) is 

\Mso+6s)) ^ -t6sa\^j,„) = riil^.J. (A27) 

We are interested in time scales long compared to 6s, 
during which many atoms interact with oscillator. There- 
fore, let us consider the time interval As — N6s, with 

^ 1, but the probability for the atom to jump to the 
excited state remains small N6s'^fiso <C 1. Thus, the 
probability that more than one jump is observed during 
As is negligible. The probability that no jump occurs 
during As is determined from the norm of the unnor- 
malized state vector conditioned on measuring each of 
a sequence of N atoms in state |0), which according to 
Eq.EMlis 

|V;,y(As)) ^^IV'.o) (A28) 
^ (l-^N6s^a^aj\iJso)- (A29) 



The probability that no jumps occur is then 

P„j(As) = (Vi„,(As)|^„j(As)), (A30) 

which simplifies to 

P,y(As) ^ 1 - N6s''{a^a)so. (A31) 

The probability that only one jump is observed at some 
point during the interval As is obtained from the unnor- 
malized oscillator state vector conditioned on measuring 
the m*'* atom in the excited state 

|^™(As)) = ri^-^l^if^^lV.o), (A32) 
-i6sa\Tpsg), (A33) 

where in the second line I have substituted in Eq. IA27I 
and retained terms only of order 5s. Noticeably Eg. IA33I 
is independent of m. With the aid of Eq. IA33[ we find 
that the probability that one jump occurs at some point 
during As, 

P,-(As) = (i"(As)|Vi7(As)) (A34) 

m=0 

- A(5s2(a^a)so. (A35) 

We now see from Eq. IA35I that the probability to ob- 
serve a transition during As is of order As = N6s, which 
is an indication that the series of jumps observed during 
the evolution of the oscillator is described by a Poisson 
process. Moreover, a majority of the time no jump will 
be observed, and the change in the state vector according 
to Ea. lA29l will be small and of order As. With probabil- 
ity Pj [Eq. lA35] though the atom will jump and the state 
vector will change dramatically under the action of the 
lowering operator a [Eq. I A33| . We may formulate this 
observation mathematically by introducing a stochastic 
increment for the Poisson process, AN^ , which is typi- 
cally zero in any small time interval As = Ai/A, but with 
probability Pj is one. Specifically, the Poisson increment 
is defined by the relations 

{AN+f = AN+ (A36) 

and 

E[AN+] = gAt{a^a)t, (A37) 
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where I have replaced s = Xt and introduced the jump 
(or decay) rate g = X^St. Using AN^ , we may combine 
Eqs.|A29]and|S3as 

A\i,t,)^-^gAta^a\i^,,)+AN+{a\ijto)-\i'to)), (A38) 
which upon taking the infinitesimal limit At dt reads 



dlM = -^gdt a'' a\ifjto) + dN^ {a\^pt^,) 



IV'*o))- (A39) 



To complete the derivation, we recognize that repeat- 
ing the above sequence of steps with an atom initially in 
the excited state leads to Eq. IA39I with the replacement 
-ir^ d. The stochastic Schrodigner equation then fol- 
lows by combining Eq. IA39I and Eq. IA39I with a' O a, 
each weighted by the likelihood that an atom is initially 
in its ground state tq or an atom is initially in its excited 
state ri , receptively. Equation [T^] is finally recovered by 
leaving the adiabatic-interaction picture through invert- 
ing Eqs. IA6I and IA17[ and normalizing the state vector 
over the small interval dt. 



Appendix B: Rotating Wave Approximation and the 
Time-Dependent Jaynes-Cummings Interaction 
Hamiltonian 

In this appendix, I argue that the interaction Hamilto- 
nian Hj{ft) in Eq. iniis a physically relevant interaction. 
I will demonstrate that Hj{ft) is the secular (or rotating 
wave) approximation of a more general interaction, and 
well approximates the dynamics on relevant time scales. 

An experimental realization of a quantum harmonic 
oscillator is an electromagnetic field confined to a su- 
perconducting microwave cavity [60l - l62j . In such exper- 
iments, the force driving the harmonic oscillator corre- 
sponds to time-dependent classical macroscopic charges 
moving determinisiticly. The two-level systems corre- 
spond to atoms traversing the cavity. Since the atoms 
are small compared to the wave-length of the electro- 
magnetic field, their interaction is well described in the 
dipole approximation. Using the notation for the har- 
monic oscillator, the atom-field interaction in the dipole 
approximation in the presence of a classical external field 
is Ml 



(Bl) 



where is the coupling strength, which depends on the 
dipole moment of the atom. Substituting in the defini- 
tions of at and a| [Eqs. [5], we find 



V{ft)^X{al+at){<J^ +a), 



(B2) 



where A = Q/^/2rnuj. 

By applying a perturbative analysis to the equation of 
motion for the coupled oscillator and atom, I will verify 



that when the coupling is weak. 



A 



e = - < 1, 



(B3) 



we can ignore the terms in V [Eq. IB2| that do not con- 
serve energy (attr and alcr'^), and approximate the evolu- 
tion by Hj{ft), at the expense of replacing at and a| by 
the displaced operators dt and aj 

The dynamical evoution of the time-evolution opera- 
tor U{t) for the coupled oscillator-atom system in the 
dipole approximation is generated by the sum of H{ft) 
[Eq.E], Ha [Eq.E] and V{ft) [Eq.EI] according to the 
Schrodinger equation 



dtU{t) = -I [H{ft) +Ha + V{ft)] U{t), 



(B4) 



with initial condition U{0) = I. F or cl arity, we switch 
to the adiabatic-interaction picture |l02l |. To obtain the 
adiabatic-interaction time-evolution operator, we apply 
the operators K\t,0) [Eg.KTO] and A'' (ft) [Eq. US] to 
U as 



UAi{t) ^ KHt,0)A\ft)U{t). 



(B5) 



A differential equation for Uai is obtained by differenti- 
ating Eg. IBSI with time; substituting in the definitions of 
K [Eq. lAlO] and A [Eq. IA3| ; exploiting their properties 
in Eqs. IA4I and IA16I and finally scaling time s = uit to 
make the e dependence explicit: 



dsUAiis) = -ieh{fs)UAi{s), 



where 



Hfs) =(a^ + a*sfl)o- + (a + as^o)a^ 

+ {a + a^,o)o-e"^*" + (a^ + a*,o)c^■'■e^'^ 



(B6) 



(B7) 



and a,,o = i(i^/a;)(e*'^^ - 1) [Eq.EIl]. 

To solve Eq. IB6I perturbatively using a two-time scale 
analysis, I introduce a slow time t = et and develop an 
asymptotic expansion for Uai in s and r. 



Uai{s) - u°{s, r) + eu^{s, t) + 



(B8) 



Replacing dg dg + edr in Eq. IB6I and substituting in 
Eq. IB8[ leads to the differential equation 



(B9) 



{ds+edr)[u''{s.T)+eu\s,T) + 



ieh{fg)[u^{s,T)+eu^{s,T) + ---] 



which we solve order by order in e in order to calculate 
the terms in the asymptotic expansion of Uai in Eq. IB8I 
Equating terms of order in Eq. IB9i we deduce the 
differential equation for , 



whose solution is 



(BIO) 



(Bll) 
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where F(t) is an unknown function of r only, with ini- 
tial value F{0) = I chosen to satisfy the initial condition 
Uai{0) = I- We fix F{t) by demanding that eu^ re- 
main smaller than vP on times s — 0{e~^) [At = 0(1)], 
maintaining the validity of the asymptotic expansion in 
Eg. IBSI up to times s — 0{e^^). To this end, we examine 
the term of order e in Eq. IB9I 

dsu\s, t) + dru\s, t) = -th{fs)u"{s, t) . (B12) 

Its formal solution, after substituting in Eqs. IB71 IB11[ 
and rn2] is 

U\S, t) =G(t) - S [z(aV + dC7'')F{T) + drF{T)] 

- i- [{e-'' - l)a + (e*" - l)cr^] F(t) 

-if dx[{a + a.^^o)ae^^''-^ + {a^ + a* o)cr'''e^'^] 
Jo 

(B13) 

where G{t) is a function of t alone with initial condition 
G(0) — 0. Clearly, at times s = 0(e^^), eu^ will be the 



same order of u", unless we set the term growing lineraly 
in s to zero: 

drF{T) = -i(aV + aa'')F{T). (B14) 

This is a differential equation for F with initial condition 
F{0) — I, whose solution is 

i^(T) = 6^*^^'**'"+"''*'. (B15) 

To lowest order in e, the approximate solution for Uai 
is obtained by substituting Eq. IB15I into Eq. IB8I and re- 
placing s = cot, 

UAi{t) ^ e"'^"^'''' \ (B16) 

valid up to times 

F(t), es = \t = 0(1). (B17) 

An identical expression to Eq. IB16l for Uai would be ob- 
tained starting from Hj{ft). Therefore, the full evolution 
of U can be approximated up to times Xt = 0(1) using 
Hi{ft) as long as X/uj <C 1. 



H. B. Callen, Thermodynamics and an introduction to 
thermo statistics, 2nd edition (John Wiley and Sons, 
New York, 1985). 
U. Seifert 
T. Schmiedl 



Phys. Rev. Lett., 95, 040602 (2005) 



T. Speck, and U. Seifert, J. Stat. Phys., 
128, 77 (2007). 

T. Schmiedl and U. Seifert, J. Chem. Phys., 126, 044101 
(2007). 

M. Esposito, U. Harbola, and S. Mukamel, Phys. Rev. 
E, 76, 031132 (2007). 

U. Seifert, Eur. Phys. J. B, 64, 423 (2008). 

M. Esposito and C. Van den Broeck, Phys. Rev. Lett., 

104, 090601 (2010). 

K. Sekimoto, Prog. Theor. Phys. Supp., 130, 17 (1998). 
K. Sekimoto, Phys. Rev. E, 76, 060103(R) (2007). 
K. Sekimoto, Stochastic Energetics, Lect. Notes Phys., 
Vol. 799 (Springer, Berlin Heidelberg, 2010). 
C. Tietz, S. Schuler, T. Speck, U. Seifert, 
J. Wrachtrup, Phys. Rev. Lett., 97, 050602 (2006). 
V. Blickle, T. Speck, L. Helden, U. Seifert, 



and 



and 



C. Bechinger, |Phys. Rev. Lett., 96, 070603 (2006) 
T. Speck, V. Blickle, C. Bechinger, and U. Seifert, Eu- 
rophys. Lett., 79, 30002 (2007). 

R. J. Harris and G. M. Schiitz, J. Stat. Mech.: Theor. 
Exp., P07020 (2007). 

C. Jarzynski, Eur. Phys. J. B, 64, 331 (2007). 

E. M. Sevick, R. Prabhakar, S. R. Wilhams, and D. J. 

Searles, Ann. Rev. Phys. Chem., 59, 603 (2008). 

C. Jarzynski, Ann. Rev. Condens. Matter Phys., 2, 329 

(2011). 

M. Esposito, U. Harbola, and S. Mukamel, Rev. Mod. 
Phys., 81, 1665 (2009). 
M. Campisi, P. Hanggi, 
Phys., 83, 771 (2011). 



[20] G. N. Bochkov and Y. E. Kuzovlev, Zh. Eksp. Teor. 
Fiz., 72, 238 (1977). 



[21] C. Jarzynski, [Phys. Rev. Lett., 78, 2690 (1997 ) 
[22] G. E. Crooks, T~Stat. Phys., 90, 1481 (19M! ' 



[23] G. E. Crooks, phys Rev. E, 61, 2361 (2000) 
[24] R. Kawai, J. M~n~ Parrondo, and C. Van den Broeck, 
Phys. Rev. L ett., 98, 080602 (2007) 



[25] C. Jarzynski, Phys. Rev. E, 73, 046105 (2006) 



and P. Talkner, Rev. Mod. 



[26] A. Gomez-Marin, J. M. R. Parrondo, and C. Van den 
Broeck, Europhys. Lett., 82, 50002 (2008). 

[27] S. Vaikuntanathan and C. Jarzynski, Europhys. Lett., 
87, 60005 (2009). 

[28] This terminology differs from that commonly used in 
thermodynamics, for example in Ref. [3], where a 
closed system exchanges energy with its environment, 
but not matter. The definition of closed (and open) sys- 
tems used in the present article conforms to the stan- 
dard usage in quantum optics -59]. 

[29] H. Tasaki, "Jarzynski relations for quantum systems 
and some ap plications," (2000), cond-mat/00092 44v2. 

[30] S. Mukamel, phys. Rev. Lett., 90, 170604 (2003)| . 

[31] P. Talkner, E. Lutz, and P. Hanggi, Phys. Rev. E, 75, 
050102(R) (2007). 

[32] P. Talkner and P. Hanggi, J. Phys. A: Math. Theor., 40, 
F569 (2007). 

[33] P. Talkner, P. Hanggi, and M. Morillo, Phys. Rev. E, 

77, 051131 (2008). 
[34] M. Campisi, P. Talkner, and P. Hanggi, Phys. Rev. 

Lett., 102, 210401 (2009). 
[35] M. Campisi, P. Talkner, and P. Hanggi, Phis. Trans. 

R. Soc. A, 369, 291 (2011). 
[36] J. Teifel and G. Mahler, Phys. Rev. E, 76, 051126 

(2007). 

[37] A. Engel and R. Nolte, Europhys. Lett., 79, 10003 



16 



(2007) . 

[38] P. Talkner, P. S. Burada, and P. Hanggi, Phys. Rev. E, 

78, 011115 (2008). 
[39] M. Campisi, P. Talkner, and P. Hanggi, J. Phys. A: 

Math. Theor., 42, 392002 (2009). 
[40] M. Campisi, P. Talkner, and P. Hanggi, Phys. Rev. E, 

83, 041114 (2011). 
[41] J. Teifel and G. Mahler, Eur. Phys. J. B, 75, 275 (2010). 
[42] G. Huber, F. Schmidt-Kaler, S. Deffner, and E. Lutz, 

Phys. Rev. Lett., 101, 070403 (2008). 
[43] G. E. Crooks, J. Stat. Mech.: Theor. Exp., 10, P10023 

(2008) . 

[44] P. Talkner, M. Campisi, and P. Hanggi, J. Stat. Mech.: 

Theor. Exp., P02025 (2009). 
[45] M. Campisi, P. Talkner, and P. Hanggi, Phys. Rev. 

Lett., 105, 140601 (2010). 
[46] F. Ritort, Physics, 2 (2009). 

[47] M. Esposito and S. Mukamel, Phys. Rev. E, 73, 046129 
(2006). 

[48 1 S. Deffner and E. Lutz, Phys. Rev. E, 77, 021128 (2008). 



[49 S. Deffner, M. Brunner, and E. Lutz, Europhys. Lett. 



[70' 
[71 

[72" 
[73; 

[74 
[75; 



[76; 

[77 

[78" 
[79^ 
[80^ 
[81 
[82" 



94, 30001 (2011). 



[83" 



[50 ] S. Deffner, O. Abah, and E. Lutz, Chem. Phys., 375, 



[51 W. De Roeck and C. Maes, Phys. Rev. E, 69, 026115 



200 (2010). 



(2004 



[84 



[52 J Y. Subagi and B. L. Hu, Phys. Rev. E, 85, 011112 
(2012). 



785 



[53 H. Carmichael, in An open quantum systems approach to 
quantum optics, Vol. 18 (Springer-'Verlag, Berlin, 1993). 
[54] H. M. Wiseman, Quantum trajectories and feedback, 
Ph.D. thesis. University of Queensland (1994). 



[86" 



"[87 

[88" 
[89" 



[55 ] M. Plenio and P. Knight, Rev. Mod. Phys., 70, 101 
(1998). 



[56 K. Jacobs, Topics in quantum measurement and quan- 
tum noise, Ph.D. thesis. University of London (1998), 
arxiv:9810015v2. 
[57] K. Jacobs and D. Steck, Contemp. Phys., 47, 279 
(2006). 

T. A. Brun, Am. J. Phys., 70, 719 (2002). 
H.-P. Breuer and F. Petruccione, The theory of open 
quantum systems (Claredon Press, Oxford, 2009). 
J. M. Raimond, M. Brune, and S. Haroche, Rev. Mod. 
Phys., 73, 565 (2001). 

H. Walther, B. T. H. "Varocoe, B.-G. Englert, and 
T. Becker, Rep. Prog. Phys., 69, 1325 (2006). 
S. Gleyzes, S. Kuhr, C. Guerlin, J. Bernu, S. Deleghse, 
U. B. Hoff, M. Brune, J.-M. Raimond, and S. Haroche, 
Nature, 446, 05589 (2007). 

T. Kist, M. Orszag, T. A. Brun, and L. Davidovich, J. 
Opt. B: Quantum Semiclass. Opt., 1, 251 (1999). 
M. Santos and A. Caravalho, Europhys. Lett., 94, 64003 
(2011). 

G. E. Crooks, Phys. Rev. A, 77, 034101 (2008). 
"W. De Roeck and C. Maes, Rep. Math. Phys., 18, 619 
(2006). 

J. Dalibard, Y. Castin, and K. Molmer, Phys. Rev. 
Lett., 68, 580 (1992). 

M. Ueda, N. Imoto, H. Nagaoka, and T. Ogawa, Phys. 
Rev. A, 46, 2859 (1992). 

S. Attal and Y. Pautrat, Ann. Inst. Henri Poincare, A, 



"[90 



[58" 
[59^ 

[60; 

[61 

[62; 

[63; 

[64 

[65 
[66; 



[67; 
[68; 
[69; 



~W 

[92; 

[93 
[94 
[95; 

[96 
[97 

[98; 

[99; 
[100; 
[101 
[102; 
[103; 



7, 59 (2006). 

S. Attal and A. Joye, J. Stat. Phys., 127, 1241 (2007). 
K. R. Parsatharathy, An introduction to quantum 
stochastic calculus (Birkhauser Verlag Basel, 1992). 
C. Pellegrini, Stoch. Proc. AppL, 120, 1722 (2010). 
S. Attal and C. Pellegrini, Open Syst. Inf. Dyn., 17, 389 
(2010). 

A. Bassi and G. Ghirardi, Phys. Reps., 379, 257 (2003). 
C. Cohen- Tannoudji, J. Dupont-Roc, and G. Grynberg, 
Photons and atoms: Introduction to quantum electrody- 
namics (Wiley-Interscience, New York, 1989). 
S. Pielawa, L. Davidovich, D. Vitali, and G. Morigi, 
Phys. Rev. A, 81, 043802 (2010). 

J. Sakurai, Modem Quantum Mechanics, edited by S. F. 

Tuan (Addison- Wesley, New York, 1994). 

T. D. Kieu, Phys. Rev. Lett., 93, 140403 (2004). 

U. Seifert, private communication. 

T. Monnai, |Phys. Rev. E, 72, 027102 (2005) . 
H.-P. Breuer, Phys. Rev. A, 68, 0321 05 (2003). 



G. N. Bochkov and Y. E. Kuzovlev, Physica, 106, 443 



(1981) 



G. Gallavotti and E. G. D. Cohen, 



2694 (1995) 



Phys. Rev. Lett., 74 



D. J. Evans and D. J. Searles, Adv. Phys., 51, 1529 



(2002) 



J. L. L ebowitz and H. Spohn, J. Stat. Phys., 95, 333 



(1999) 



T. Hat ano and S. I. Sasa, Phys. Rev. Lett., 86, 3463 



(2001) 



C. Maes, Seminaire Poincare, 2, 29 (2003). 

R. D. Astumian, Phys. Rev. E, 76, 020 102 (2007). 



C. Maes, F. Red ig, and A. V. Moffaert, [J. Math. Phys 
41, 1528 (2000)1 . 

C. Maes and K. Netocny, |j. Stat. Phys., 110,^9 



(2003) 



T. M. Cover and J. A. Thomas, Elements of Information 

Theory, 2nd ed. (Wiley-Interscience, 2006). 

I. Callens, W. De Roeck, T. Jacobs, C. Maes, and 

K. Netocny, Physica D, 187, 383 (2004). 

H. M. Wisema n. Quantum Semicl. Opt.. 8. 2 05 (1996). 

G. E. Crooks, |Phys. Rev. E, 60, 2721 (1999)| . 

R. B. GrifSths, Consistent Quantum Theory (Cam- 
bridge University Press, United Kingdom, 2002). 
T. A. Brun, Phys. Rev. A, 61, 042107 (2000). 

H. Maassen, in Quantum Probability Communications 
XII, edited by S. Attal and J. M. Lindsay (World Sci- 
entific, Singapore, 2003) pp. 23-58. 

H.-S. Goan, G. J. Milburn, H. M. Wiseman, and H. B. 
Sun, Phys. Rev. B, 63, 125326 (2001). 
K. Jacobs, P. Lougovski, and M. Blencowe, Phys. Rev. 
Lett., 98, 147201 (2007). 

S. Peil and G. Gabrielse, Phys. Rev. Lett., 83, 1287 
(1999). 

D. Spehner and M. Orszag, J. Math. Phys., 43, 3511 
(2002). 

S. Klarsfeld and J. A. Oteo, Phys. Rev. A, 45, 3329 
(1992). 

P. Pechukas and J. C. Light, J. Chem. Phys., 44, 3897 
(1966). 



